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We present an efficient method to solve scattering problems in two-dimensional open billiards with 
two leads and a complicated scattering region. The basic idea is to transform the scattering region to 
a rectangle, which will lead to complicated dynamics in the interior, but simple boundary conditions. 
The method can be specialized to closed billiards, and it allows the treatment of interacting particles 
in the billiard. We apply this method to quantum echoes measured recently in a microwave cavity, 
and indicate, how it can be used for interacting particles. 
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I. INTRODUCTION 

The role of two-dimensional billiards in the develop- 
ment of chaos theory in general and of wave chaos in 
particular is basic. The pioneering work of Sinai Q has 
opened the door for a large field of research where work 
of Bunimovich and others H has provided us with 
a number of interesting shapes, for which chaoticity (at 
least ergodicity and the Kolmogorov property) have been 
proven. This has become particularly relevant as for soft 
potentials, at this point, no such example was found. 
There is even a wide-spread belief, that no example may 
exist. The field got additional impulse, from the fact that 
in scattering systems hypcrbolicity may be easily proven 
e.g. for the two and three disk problem p|, and for the 
latter there exists a bounded manifold, that is chaotic. 
For scattering systems this property is also achievable 
for soft potentials [f|. The study of billiards has also 
proven interesting in mixed systems, both of the bound 
(e.g. H0) and open type H|^. 

In physics, billiards have acquired increasing impor- 
tance as they are seen to represent relevant aspects of 
systems as small as quantum dots E3 and as large as 
planetary rings E3- Wave realizations of such objects 
have become even more popular, as experiments with flat 
microwave cavities become available from an ever larger 
number of labs H3, E3, E3, EH EH E3, EM EJ] - Such 

experiments emulate single particle properties of meso- 
scopic devices. Furthermore such billiards provide a very 
direct and experimentally accessible way to test hypoth- 
esis on classical-quantum connections, such as the rela- 
tion between the classical ensembles of random matrices 
[13 and chaotic evolution, for spectra [M0E1EHE1, 
wave- functions [HE3 and scattering properties |3fll3l| . 

Microwave billiards emulate the wave behavior of the 
corresponding classical billiard and thus they form the 
most direct way to test hypotheses proposed for classical- 
quantum or ray-wave connections. For closed billiards 
extremely efficient algorithms are exist [31 [33| , while for 
scattering situations only some standard finite element 
programs are available to solve general problems. Re- 
cently a new technique has been proposed for special sit- 
uations, such as the ripple channel |38ll39|. The purpose 



of the present paper is to generalize this technique such 
that it can bee applied to typical open billiards in use as 
quantum dots or in microwave experiments. The basic 
idea is to transform the complicated border to a simple 
rectangle with openings. The price we pay is that we now 
have to deal with a much more complicated Hamiltonian 
instead of the simple Laplace-Beltrami operator. 

The advantage of this procedure thus clearly results 
from the simple boundary conditions. While these also 
play a role for bound systems they are more pronounced 
in scattering systems. If we use R-matrix theory we have 
to implement Dirichlet and Neumann conditions in differ- 
ent parts of the surface simultaneously, and this becomes 
rather easy with the simple boundary conditions. The 
more complicated dynamics can be handled by in such 
a way, that the main numerical effort is reduced to fast 
Fourier transforms, which are known to be extremely ef- 
ficient. 

We apply our method to obtain transmission and the 
power spectrum for a scattering configuration with mixed 
hilbertspace in the interior of the cavity very similar to 
the one for which scattering experiments have been per- 
formed successfully in Darmstadt with a superconducting 
microwave cavity |8j . We shall refer to this configuration 
as Darmstadt experiment configuration (or billiard) and 
we will show that we can reproduce the main feature, 
namely the scattering echo [ij, which was theoretically 
predicted. Also we find, that this kind of dynamics intro- 
duces a step-structure in total transmission, which should 
be observable in microwave billiards, though total cross- 
section measurements over many channels are always dif- 
ficult. Furthermore we show, that two-body interactions 
can be introduced efficiently, thus opening the door for a 
plethora of additional important applications. 

We shall proceed in the next section to recall some el- 
ements of S- and R-matrix theory. Furthermore we shall 
connect these to concepts such as reflection and transmis- 
sion, much more common in the language of mesoscopics. 
In the section III we attack the main problem, namely to 
show how we transform a complicated shape to a rectan- 
gle and what operator will appear in calculation of the 
the R-matrix. We then obtain a general result with a 
well defined numerical recipe, which is mainly based on 
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fast Fourier transforms (FFT). In section four we apply 
this technique to the Darmstadt experiment billiard used 
in [3] to experimentally confirm the existence of scatter- 
ing echoes due to large integrable islands in the classi- 
cal system. Note that the original theoretical prediction 
for a quantum system was made generically using a sim- 
ple delta-kicked one-dimensional system, and the authors 
there claim, that a computation of the scattering billiard 
is not feasible to sufficiently short wave length with avail- 
able means 0- In the following section we open the door 
to further applications as we discuss, how the method 
can be applied for systems with two-body interactions. 
We then give some brief conclusions. 
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FIG. 1: (a) An open billiard of length L . (b) the scattering 
cavity is transformed to a rectangular region in u-v plane 



II. SCATTERING IN OPEN BILLIARDS 

The system we are concerned with is sketched in Fig. 
n It consist of a flat wave guide of width w that sup- 
ports iV open modes (or channels) with a rather general 
scattering region of length L. The specific method we de- 
velop, requires two leads (openings are in principle also 
possible) at opposing ends of a scattering region. We 
shall arbitrarily fix them to be on the right and on the 
left, as in the figure. The scattering region between n the 
leads or openings shall be described by two single- valued 
functions of a coordinate say x defined between the two 
leads. No additional scatterers in the so defined cavity 
are allowed. We will assume, that the leads have equal 
width. The latter is not essential for the method but sim- 
plifies notation. The scattering problem will be defined 
in terms of a 2N x 2N scattering matrix S. Having nano 
structures in mind we can rewrite the S-matrix as 



S 



t 



(1) 



where r (r') and t (t 1 ) are the reflection and transmission 
matrices for incidence on the left (right) of the scattering 
region. The dimensionless conductance is obtained from 
S as 



T = tr(tt f ) 



(2) 



Our calculations will be based on the traditional R- 
matrix approach |34l l35| . This formalism has been 
adapted to obtain the S matrix for a wave guide in previ- 
ous works j3|| l38l l39| . Until now only a restricted choices 
of scattering regions has been investigated. For example 
only one of the wall was taken of sinusoidal shape while 
the other one was kept flat (ripple billiard). We want to 
generalize this approach to a generic open two-lead cav- 
ity without obstacles, where one can choose any unique 
differentiable function of P(x) and Q{x) in Fig. 1 to rep- 
resent the upper and lower walls of the scattering cavity. 

We wish to remind the reader, that the i?-matrix ap- 
proach relates the S-matrix to a matrix R determined by 



boundary conditions relating internal and external func- 
tions at the edges, where the leads are connected to the 
scattering region, such that 



^|r = i?V^| r 



(3) 



at x — r as in of the boundaries represented as dotted 
lines in Fig. 1. i?-matrix can be generalized for number 
of lead connection and number of modes in the leads. 
i?-matrix can be related to .S'-matrix, 



S 



1-ikR 
1 + ikR 



(4) 



Rather than recalling details of i?-matrix theory we 
will evoke the spirit of the method by presenting a trivial 
one dimensional example: 

In Fig. 2 we show the ID example we want to solve us- 
ing i?-matrix method. In this example we have 2 asymp- 
totes as in the general case we discuss but the r and t 
are functions instead of a matrices as in the 1 mode case 
for the real cavity problems. The exact S'-matrix can 
be found in this case and the transmission probability is 
given by 



1 



(l + ( 



(E(E-Va))(2sm(k))' 







(5) 



for a constant potential step of height Vq and length 1. 

Numerically we use a basis for the Reaction region, (R 
in Fig. 2), which is given by cosimirx), m = 0, 1, . . . oo. 
xi = and x r = 1 have been chosen. Using this basis, 
the R-matrix elements are given by, 



Rrl = 



1 oo 



E-V ' ^ E - m 2 7r 2 - V 

m—l 

2 cos(rn7r) 
E - Vq ' ^ E - m%2 - v Q 

rn—1 



E 



Ru 



= Rlr- 



(6) 



This series should be truncated at some finite value for 
a numerical calculation; m — 1000 is used in Fig. 2. The 
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FIG. 2: Scattering from a ID barrier. The solid line shows 
the transmission probability as a function of energy and stars 
are the numerical calculations. In the inset the geometry of 
the scattering system is shown. / indicates the interaction 
region, A the asymptotic regions. 



S matrix is obtained from R matrix as 



conditions at the open boundaries, and Dirichlet condi- 
tions at the walls of the billiard. With other words, the 
derivative of the wave function is set to be zero at the 
dotted line in Fig. 1 and the wave function itself is set to 
zero at surface of the billiard. Thus we have the equation 



H\Hf(x,y) >=-^l- 



■V 2 y )^(x,y) = E n ^(x,y), (7) 



with the boundary conditions V x ^(a;, y)\x=o 



\x=L 



0, ^{x,y)\ y=P{x) = 



= 0, 
and 



i &(x,y)\y=Q(x) = 0. The eigenvalues and eigenvectors 
of this equation are used to represent the scattering dy- 
namics inside the cavity. 

The solution of the Schrodinger equation in the coor- 
dinates x,y is complicated due to the arbitrary shape of 
the boundaries so we use following transformation to new 
coordinates u, v. 



S = 



1 

t 





— ik 



1 o 
o 1 



ik 



Rlr 



Rrl 

Ru 



1 
1 



ik 



Rrr Rrl 
Rlr Rll 



1 

e 



- i k 



Using this S matrix we plot the transmission probability 
,|5i,2| 2 hi Fig. 1. We note that generalizing the potential 
in this example to any shape can be done by using a finite 
difference method. |40| 

The method to obtain the S-Matrix thus proceeds in 
two steps: First obtain the solution of the boundary value 
problem inside the scattering region with Neumann and 
Dirichlet boundary conditions at lead connections and 
cavity walls respectively. Second connect the known so- 
lution in the asymptotic region(leads) to the one obtained 
for the scattering region by imposing continuity of the so- 
lution. In the method we present surprisingly the most 
time consuming step in numerics will be to obtain a com- 
plete set of function to expand the scattering wave func- 
tion inside the cavity region, which is a prerequisite to 
the first step. 



III. THE METHOD 

We shall thus start by presenting a very efficient 
method to deal with the problem inside the scattering re- 
gion. This method can also be applied to closed cavities, 
as we will point out further down. As we will use trans- 
formed coordinates the second step will not be entirely 
trivial and will be presented once we have the internal 
solutions. 



y ~ Q(x) 
P(x) - Q{x) 



(8) 



This change of variable transforms the complicated 
boundary to a simple one, i.e. a rectangular region. The 
price we pay, is that the form of the Schroedinger equa- 
tion becomes more complicated. The Hamiltonian oper- 
ator in the u—v plane can be obtained by using the above 
transformation and complete derivatives. The final result 
is 



H = Vf, 



J 



V v + 2^ iQu+ j V 2 Ju)Ju ) V..(?)i 



where J is defined as J — P — Q and subscripts indicate 
partial derivatives. In this form the symmetric nature of 
Hamiltonian is not obvious. Therefore we will represent 
the same equation in the following matrix form,[4l| 

H=j(y a Jg a ' p V(i) a = u,v (3 = u,v (10) 
where g is a metric and equivalent to, 



J 



-(Q u +vJ u ) 
J 

i+(Qu+vJ u ) 2 
7 1 



(11) 



The meaning of J becomes clear when we notice it is 
related to the determinant of the metric. It plays the 
role of Jacobian of the transformation from the x-y to 
the u-v plane. The wave function in the u-v plane can 
be expanded in terms of a basis 



A. Cavity Region 



We want to solve the Schrodinger equation in the cav- 
ity with specific boundary conditions, namely Neumann 



Eh 12 ///tvru, . .. . 

"' m 7j 7l COS(_ zT ) sin(n7rv)(12) 
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which satisfies the boundary conditions at u = 0, u — 
L and v = 0, v = 1 automatically. To determine the 
unknown coefficients B\ m we will use the orthanormality 
condition. In the original coordinates plane waves are 
orthonormal; in terms of u, v they must be orthonormal 
with a weight function given by the metric 



dudvJ^k(u,v)^k'(u,v) = Sk,k> (13) 



After multiplying the eigenvalue equation with another 
eigenfunction and integrating we obtain the matrix equa- 
tion, 



Bf 



(14) 



where we use new indices such that (n, m) — > I and 
(n',m') — > V . Eigenvalues E k correspond to energy and 
eigenvectors contain the unknown coefficients B\ of the 
wavefunction. Using the matrix form of H we have, 



1 , 



(15) 



We can integrate by parts. The surface terms are zero 
due to the boundary condition we have chosen. We thus 
find 



WMg^Vp^ = E k S, 



(16) 



Since the metric is symmetric, the matrix obtained from 
this equation is real symmetric and positive definite 
which guarantees positive energy eigenvalues. 

Once the eigenvalue equation is solved one can obtain 
scattering wave function for any energy within the range 
of wave length we allow. This makes the method efficient 
compared to direct integration of the differential equation 
by some mesh based method like finite element j3g] ■ The 
cut off on the number of basis states determines the accu- 
racy and range of validity of the method. In practice it is 
necessary to decide how many modes to choose in u and v 
direction. A corresponding integrable geometry gives an 
idea about these numbers. Also one needs to choose only 
the lower part of the eigenvalue spectrum in the bound 
region to ensure that truncation wont bring errors due 
to poor solutions of the bounded problem. For example 
m = 90 , n = 50 and using the first 3000 eigenvalues 
of the 4500X4500 matrix is suitable for a rectangular re- 
gion with length twice its width. The limit of the highest 
wavelength that can be attained in leads is connected to 
the lead size to the scattering size. If the size of scatter- 
ing region is much bigger than the size of the leads than 
the method is limited to only couple of modes in leads 
due to the cut off on the number of eigenvalues that has 
to be done in cavity region. In the application we discuss 
in this paper this is not the case since size of cavity is 
smaller than leads. 

We need an efficient method to calculate the matrix 
elements. Each matrix elements is a double integral over 



u and v. Due to the form of the metric and sinusoidal 
nature of the basis states we can achieve this writing the 
terms of Eq. ^J|as follows: 



< 9\H\Hf>= j duW u &) -7V«(%) • J dvS n S n , 

+ n'n J dttV u (-^)(— Q„)-^= ■ J dvS n C n > 

+ n'n j ' duW u (^^)(-J u )-jj ■ j dvS n vC n > 

+ mr J du-^(—Q u )-j= ■ J dvC n S n > 

+ mi J du—j=(—J u ) ^J— ■ J dvC n vS„i 

+ nnV 2 J duC m C m , (1 + 2 Qtt) ■ ! dvC n C n , 

+ nnn 2 [ duC m C m i - — • / dvC n vC n i 



J 2 



+ nn! ■k' 1 I duC m C m i 



j 2 



dvC n v 2 C n >. 



Here we use the abbreviations C m — > cos(rn7ru / L) , S n — > 
sin(rt7rw) and C n —> cos(nirv). The integral over v can 
be performed analytically. All the terms except the first 
one are of the form of a function multiplying a C m term. 
These integrals can be done using a fast Fourier trans- 
form. The first term can be put in the same form, after 
some algebra: 



duC m \7 u {C' m )V u {\og(J)/2) 
duC m .W u {C m )V a {\og{J)/2) 
duC m C rn ,(V u (log{J)/2)) 2 ) 



Therefore, when the functions P, Q and J are given as 
an array of dimension M, the integral is equal to 



J duC m F{u) = L Re ^ FFT {F) m _ 1 ) (17) 



i.e. the m — 1th element of the real part of the Fourier 
transform of the function is equal to the integral of the 
same function multiplied by C m . Note that one needs to 
incorporate end point correction [4^ | to be more accurate; 
one way is to extend integration function symmetrically 
at the end point and halve the result. This is very effi- 
cient since each array of integrals is calculated by a single 
Fourier transform. A total of seven Fourier transform is 
enough to calculate all integrals for this problem. The 
overall matrix elements can be written using the result 
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of exact integration in v direction as 

.(TO - 1) 2 7T 2 



< ¥|JI|tf >=<W(- 



L 2 



—Fi(m, to) + F2{m, to') 
+Di(n', n)F^,{m! , to) + D2(n, n')F^{m l , to 
+L>i (n, n')F 3 (m, to') + L> 2 (", «')^4 to' 
+1)3(71, n')F§(m, to') + Di{n 1 n')FQ(m, to' 
+£> 5 (n,n')F 7 (TO,TO') 



where i> integrals are obtained analytically, 

(-1 + (-!)("+"') ' 



(18) 



Di(n 7 n') = (1 - 5 n:n i)n'n 



D 2 (n, n) 



+ (1 - K,n>)nn' 



(n 12 — n 2 ) 

_^(n+n') 



(n 12 — n 2 ) J 

D 3 (n,n') = 6 n , n ,(n\ 2 /2) 
D 4 (n,n') = <5„,„'(nV/4) + 
, /(-l + (-l)(""-"))(n 2 +n' 2 ) 



D 5 (n, ri) = <5„, n ,(nV/2)(l/6 + l/(4nV)) + 
, /2(-l)( rm -")(n 2 +n' 2 ) 



(1 ~ &n,n')nn 



(ri' 2 — n 2 )) 2 



(19) 



and u integrals 



Fx (to, to') = -j^Re(FFT m (y u {C m ,)V u (log(J/2)))) 
F 2 {m,m') = ^Re{FFT m {C m ,Wu{log{J/2)))) 
F 3 (to,to') = ^Re(FFT m (-Q u V u (^))) 

F 4 (to,to') = ^Re{FFT m {-J u V u {^))) 

F 5 (m,m') = ±Re(FFT m ( iCm ' + £ m ' Ql) )) 
F 6 (to,to') = J±Re(FFT m (2C m ,Q^)) 

F 7 (m,m') = -L Re (FFT m (C ml ^)) 



(20) 



B. Leads and Scattering Matrix 

With the knowledge of states in the cavity region it 
is possible to couple known asymptotic, such as leads to 
the cavity. The corresponding solution on the left and 
on the right leads to 



l n ik n x 



1 ^sin(^) 
w 

nity > 



n e lk " x I sin 



w 



(21) 



where w is the lead width, a n , b n ,c n and, d n are scattering 
amplitudes and the wave vector is given by 



/nir\ 2 2m_E 



V w / 



h 2 



(22) 



We here assume propagating modes but it is possible to 
add evanescent modes with complex k vectors to the cal- 
culation. Their effect is discussed in ref. [3jj 

As shown in references the energy cigcnfun- 

tions \E > of the total scattering system have contribu- 
tions from both cavity and leads. 



< x,y\E >= V -yjM*>v) + E( r «*2 + r W) (23) 



3=1 



n=l 



where <j>j are basis states calculated for the cavity alone. 

The continuity of scattering wave functions at the lead 
boundary gives the connection between the leads and the 
cavity region: 



00 

E 

n' = l 
00 

£ 

n' = l 



R aL (n,n')V x *«\ XL 
R aR {n, ri)V x ^\ XR 



(24) 



where 



R a p{n,ri) = —J2 



H 2 ^2, <f>j tn (x a )<l>j in >(X0) 



E -is 



(25) 



The final step is to find eigenvalues and eigenvectors of 
this matrix which in turn gives yield the unknown expan- 
sion coefficients of the wave function in the u-v plane and 
obviously the energy eigenvalues. It is possible to obtain 
the wave function as functions of x, y applying the in- 
verse transform to the wave function in in terms of it, v. 
Since the expansion coefficients are already determined 
this transformation is straight forward. 



is the (n,n')th matrix element of the R matrix. Here 
4>j^ n is the overlap between cavity states and lead func- 
tions and given by 



/2 f 00 niry 
- / dy<pj (x a )s'm( ) (26) 
J-J .In W 



where a = L, R and xl = 0, Xr = L in Fig. 1. 
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FIG. 3: (a) The Darmstadt- experiment geometry of length 
L = 4.32 with leads W = 1. (b)same as in (a) but with some 
small disorder in the middle of the billiard, (c) The lower wall 
in (a) has been replaced by a disordered wall. 
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FIG. 4: Conductance versus wave vector for 3 different ge- 
ometries corresponding to (a), (b) and (c) in Fig. 3 



IV. APPLICATION: QUANTUM ECHOES 
FROM CAVITIES WITH A LARGE 
INTEGRABLE ISLAND 

As an application of our method we solve the scattering 
problem in the Darmstadt-experiment geometry shown 
in Fig. 3. A similar geometry to Fig. 3(a) has been used 
in recent experiments Q . The analytical curves of upper 
and lower wall in Fig. 3a(a) are given as follows, 

(XX^ f ^x^ \ 
P(x) = Acxp(--^-) and Q(x)=\i(3 'JY ) ■ 

We use the same parameters as in experiment such as 
a = 0.161, (3 = 0.2 and 7 = 0.1. While in experiment 
a scaling factor of A = 5cm is used we use A = 0.432 in 
arbitrary units and we shift the origin of the coordinates 
to the beginning of the left lead for convenience. Note 
though, that the original experiment does not connect 
to leads! Rather the opening is arbitrary, though in fact 
naturally finite. We will return to these differences later. 

Apart from the original geometry we also consider a ge- 
ometry with some small irregularity in the center of the 
cavity. Such an irregularity should destroy the stability 
of the central periodic orbit in the classical system and 
thus at least damage or deform the large regular island 
significantly. This corresponds to the mentions introduc- 
tion of a perturbing scatterer in the experiment, which 
destroyed the echo [8j . Furthermore we consider that the 
entire parabolic wall may perturbed in a disordered way. 

After calculating the S matrices for these systems we 
found the total conductance as described in section 1. 
The results for the 3 different geometries can be seen in 
Fig. 4. We see that for the geometry used in experiment 
there are steps in the conductance. These steps do not 
correspond to the opening of channels in the leads. This 
is because of the smaller width of the scattering region. 
The minimum width in channel is w m i n — 0.3097 with 
corresponding k e ff =3.23 matches where the first step 
seems to start in Fig. 4. With other words, this is the 
point where we pass from evanescent or tunneling trans- 
mission to open channels in the necks presented by the 



Darmstadt-experiment configuration. While the initial 
rise thus is rather trivial, the step structure is not, as 
is sensitive to even a small disturbance. The solid line 
in Fig 4 shows the total conductance through the same 
geometry but with a small wiggle at the center of the cav- 
ity. The exact shape of the wiggle is not very relevant. 
We choose a sinusoidal perturbation of the upper wall in 
the form 0.01 sin(107rx s /L) where x s is the x coordinate 
in the interval, 9L/20 < x s < 11L/20. (In Fig. 4 a 10 
times bigger scaled version of perturbation is also shown 
to give an idea of the shape of perturbation.) The effect 
of perturbation is huge on conductance data. The step 
structure tends to be lost for higher energies. 

We also show the result of surface disorder on conduc- 
tance by the dot-dashed curve in Fig. 4. We implemented 
surface disorder by first dividing the boundary curve into 
n pieces and move these pieces by a random amount, n in 
y direction. We apply a spline interpolation to connect 
the pieces to form a smooth boundary. We see that for 
n = 0.2 and n = 100 the step structure is lost and only 
resonance transmission can be seen. We note that a par- 
tial step can be observed for small disorder parameters. 

We calculate the length spectrum by taking a Fourier 
transform of the transmission amplitude 

t nm (L) = / dk F t nm {k F ) exp(~ik F L). (28) 

Jk min 

This identification can be done based on a semi-classical 
description where one can sum over different periodic or- 
bits to get the transmission amplitudes ^3| . In Fig. 5 we 
show this quantity for t\\{L) over an interval 1 < k F < 19 
with k in units of n/w. Note that in this interval we 
have a an increasing transmission matrix size each time 
a channel opens, we added the first element of each 
transmission matrix to form t\\ over the given k inter- 
val. In Fig. 5a we see the periodic oscillations starting 
at L — 4.32. Since length can be associated to time for 
a constant speed, such as expected in an electromagnetic 
cavity, this is similar to the echos observed in experi- 
ment. The period seems to be constant in the data with 
a second harmonic entering at higher lengths suggesting 
a shortening in the period. In the experiment a slight 
shortening of the period for short times appears. This we 
do not see. On the other hand the experiment does not 
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FIG. 5: (a) The length spectrum for the geometry in 
Fig. 4(a), (b) same for the geometry in Fig. 4(b) (c) same 
for the geometry in Fig. 4(c). 



show the drastic shortening with the higher harmonics. 
We mus keep in mind, that the problems are not iden- 
tical, because of the finite diameter leads we attached 
here. Also we measure transmission to the entire channel, 
rather than to an antenna placed near a semi-cclassically 
significant manifold as the latter is not meaningful in a 
billiard with leads. On the other hand effects of trun- 
cation must be studied carefully in each case. At this 
point we wish to stress, that the semi-classical structure 
associated with the scattering echoes again is easily seen, 
and - as in the experiment - for wave length far from 
semi-classics. 

To get a better impression we look at the power spec- 
trum for the geometry in Fig. 3(a). In Fig. 6(a) we show 
this spectrum calculated from k(-K/w) = 6 to k(ir/w) = 9. 
Since the minimum transmission matrix is 6X6 we sum 
over 6 modes. This region roughly corresponds to the 
third plateau in Fig. 4. 



(29) 



where t nm is calculated using equation 1281 but with ap- 
propriate limits for k integration. In this figure we ob- 
tain peak values at L\ = 4.78, L 2 — 6.65, L3 = 8.05 
,L 4 = 9.1,^5 = 10.2. The peak value corresponds to the 
simplest direct transmission orbit starting at the middle 
of the channel opening and hitting the center of Gaussian 
shaped wall before exit. As it can be seen here period 
decreases with increasing length. We observe oscillations 
up to high lengths as can be seen in inset. 

We also show similar calculations in the interval 
k(n/w) — 10 < k(n/w) < k(n/w) — 13 including 10 
modes in Fig. 6(b). Although we have different limits to 
the integration and different size of the matrix. We get 
results very similar to those the Fig. 6(a). 
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FIG. 6: (a) Power spectrum from k(ir/d) = 6 to k(w/d) — 9 
including 6 modes. (b)Power spectrum from k(n/d) = 10 to 
k(ir/d) = 13 including 10 modes. 



V. TOWARD MANY-BODY PROBLEMS 

Real naosystems have interacting particles, and treat- 
ment of such systems is incipient. We shall show, that 
the methods presented here is very promising for such 
situations, at least as long as the particles enter or leave 
the quantum dot only one at a time. The basic problem 
he consists of treating two-body interactions, as we may 
typically encounter in a few-body system. The two par- 
ticle inside a wave guide will provide the main difficulty. 
It is also an important system on its own right (see the 
review in |44|). 

Lets consider two particles inside a cavity region as 
shown in Fig. 1. We will omit the problems asso- 
ciated with anti-symmetrization or symmetrization for 
Fermions or Bosons respectively, by assuming distin- 
guishable particles, but standard methods can be used 
to take care of the corresponding symmetry adaption as 
necessary. The total Hamiltonian is, 



Vl 



-V 



y-2 



■V(\xi-x 2 \,\vi-V2\) (30) 



with a potential V which depends on the distance of the 
two particles. The interaction term can be written in 
terms of single particle wave functions as follows, 



Hijkl — 
V (I 



J J dxidx 2 dyidy 2 (j) l (xi,yi)4 l j (x 2 , y 2 ) 

V2\)(f)k(xi,yi)(l>i(x 2 ,y 2 ) (31) 



xi - x 2 \, \m 



where (f>i(x,y) is the single particle solution in the cav- 
ity. Energies of the interacting particle can be found after 
diagonalizing this matrix. As the form of the single parti- 
cle functions is not simple this is not an easy calculation. 
But as we discussed in this paper, we can transform this 
integration to u-v plane where we know the expansion 
of the single particle eigenfunctions in terms of simple 
harmonics. In u-v coordinates the distance between par- 
ticles wont be a line but a curve given by the metric. The 
matrix element in u-v coordinates becomes, 

Hij k i = J J duldu2dvldv2(f> i (ul,vl)(j) j (u2,v2) 

V {ul,u2,vl,v2)<i) k {ul,vl)<l) l {u2,v2) (32) 
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where <p l (u,v) = ^2 mn B l m n cos(miru/L) sm(nirv) with 
coefficients B l mn calculated from a single particle equa- 
tion. The harmonic expansion makes the integration sim- 
ple. 

VI. CONCLUSION 

We present a new method to treat open billiards, that 
has significant advantages over older techniques. It is de- 
signed for billiards with two leads, though generalizations 
are possible. The basic limitation is that both irregular 
boundaries that connect the leads be single valued func- 
tions of some space coordinate, that connects the center 
of the two leads. We then transformed the scattering re- 
gion into a rectangle and calculate the quantum dynamics 
which become involved. The boundary conditions now 
result trivial, and this allows to calculate the R-matrix 
easily and thus solve the scattering problem. It is possible 
to write the dynamics of two or more interacting particles 
in this form, as long as we do not attempts to transform 



to relative and center of mass coordinates. This would 
destroyOy the simple boundary conditions. To show the 
power of the method we calculated the quantum scatter- 
ing echoes for such a two lead structure very near to one 
actually used in a microwave experiment. The results 
of this experiment were not accessible to a numerical 
calculation with standard microwave programs, but we 
obtained results essentially consistent with experiment, 
though we see a secondary effects not apparent in the 
experiment. Future research will have to show whether 
this effect is suppressed in the experiment, whether it is 
due to the use of leads or whether it is a consequence of 
the cutoff we use. This should not distract from the fact, 
that the main effect can be seen clearly. Concerning the 
example we also note that the calculation reveals a step 
structure in the transmission, which is directly related to 
the scattering echoes and was not previously noted. 
Acknowledgments 

The authors are grateful for useful discussions to 
C. Jung. Financial support under DGAPA project 
IN191603 and CONACyT project 43375 is acknowledged. 



[1] Ya. G. Sinai, Russ. Math. Surv. 25, 137 (1970). 

[2] L.A. Bunimovich, Funct. Anal. Appl. 8, 254 (1974); 

[3] L. Bunimovich, G. Casati, and I. Guarneri Phys. Rev. 

Lett. 77, 2941 (1996) 
[4] B. Eckhardt J. Phys. A 20, 5971 (1987) 
[5] C. Jung and H.J. Scholz, J. Phys. A 20, 3607, (1987); C. 

Jung and P. Richter J. Phys. A 23, 2874 (1990) 
[6] L. A. Bunimovich, Chaos 11, 802808 (2001). 
[7] Lansel, Steven; Porter, Mason A.; and Bunimovich, 

Leonid A. Chaos, Vol. 16, No. 1 013129. (2006) 
[8] C. Dembowski, B. Dietz, T. Friedrich, H.-D. Graf, A. 

Heine, C. Mejia-Monasterio, M. Miski-Oglu, A. Richter, 

and T. H. Seligman Phys. Rev. Lett. 93(13) 134102, 

(2004) 

[9] C. Jung, C. Mejia and T.H. Seligman, Europhys. Lett. 
55 616 (2001) C. Jung et al., New J. Phys. 6, 48 (2004) 
[10] C. M. Marcus, A. J. Rimberg, R. M. Westervelt, P. F. 
Hopkins, and A. C. Gossard, Phys. Rev. Lett. 69, 506 
(1992). 

[11] L. Benet and T.H. Seligman, Phys. Lett. A 273, 331 

(2000); L. Benet, J. Broch, O. Merlo and T. H. Seligman, 

Phys. Rev. E. 71 036225 (2005). 
[12] H. J. Stockmann and J. Stein, Phys. Rev. Lett. 64, 2215 

(1990). J. Stein and H. J. Stockmann, Phys. Rev. Lett. 

68, 2867 (1992). 
[13] Y.-H. Kim, M. Barth, H. J. Stockmann, J. P. Bird, Phys. 

Rev. B 65, 165317 2002. 
[14] H. D. Graf, H. L. Harney, H. Lengeler, C. H. Lewenkopf, 

C. Rangacharyulu, A. Richter, P. Schardt, and H. A. 

Weidenmuller, Phys. Rev. Lett. 69, 1296 (1992). 
[15] C. Dembowski, B. Dietz, H. D. Graf, A. Heine, T. Pa- 

penbrock, A. Richter, and C. Richter, Phys. Rev. Lett. 89, 

064101 (2002). 

[16] R. Hofferbert, H. Alt, C. Dembowski, H. D. Graf, H. L. 
Harney, A. Heine, H. Rehfeld, and A. Richter, Phys. Rev. 
E 71, 046201 (2005). 



[17] E. Doron, U. Smilansky and A. Frenkel, Phys. Rev. Lett. 

65, 3072 (1990) 
[18] S. Deus, P.M. Koch and L. Sirko, Phys. Rev. E 52,1146 

(1995) 

[19] S. Shridar, Phys. Rev. E 67, 785 (1991) 

[20] D.H. Wu et. al. Phys. Rev. Lett. 81, 2890 (1998) 

[21] J. Barthlemy, O. legrand and M. Montessagne, Phys. 

Rev. E 71,016205 (2005) 
[22] E. Cartan, Abh. Math. Sem. Univ. Hamburg, 11 116 

(1035) 

[23] G. Casati, F. Valz-Gris and I. Guarneri, Lett. Nuovo Ci- 
mento 28, 279 (1980) 

[24] O. Bohigas, M. J. Giannoni, and C. Schmit Phys. rev. 
Lett. 52, 1 (1984) 

[25] M. V. Berry, Chaotic Behaviour of Deterministic Sys- 
tems, edited by G. Iooss, R. Helleman, and R. Stora 
(North- Holland, Amsterdam), (1983), 

[26] F. Leyvraz and T. H. Seligman, Phys. Lett.A 168, 348 
(1992) 

[27] S. Mller, S. Heusler, P. Braun, F. Haake, and A. Altland 

Phys. Rev. Lett. 93, 014103 (2004) 
[28] C.E. Porter, R.G Thomas, Phys. Rev. 104, 483 (1956) 
[29] T Gorin, Tomaz Prosen, T. H. Seligmann New J. Phys. 

6 20, (2004) 

[30] J. Verbaarschot, H.A. Weidemueller and M. Zirnbauer, 

Phys. Rep. 129 3672 (1985) 
[31] C. Jung, C. Lipp, and T. H. Seligman, Ann. Phys. (NY.) 

275, 151 (1999) 
[32] E. Vergini and M. Saraceno, Phys. Rev. E 52, 2204 

(1995). 

[33] Tomaz Prosen, Phys. Rev. E, 60 1658 (1999) 
[34] E.P Wigner and L. Eisenbud, Phys. Rev. 72, 29 (1947) 
[35] A.M. Lane and R. G. Thomas, Rev. Mod. Phys. 30, 257 
(1958) 

[36] L. E. Reichl The transition to Chos, Second Edition 
(Springer- Verlag, 2004.) 



9 



[37] H. Feshbach Topics in the Theory of Nuclear Reactions, 
in Reaction Dynamics, (Gordon and Breach, New York, 
1973). 

[38] G. B. Akguc and L. E. Reichl, Phys. Rev. E 67, 46202 
(2003). 

[39] G. B. Akguc and L. E. Reichl, Phys. Rev. E 64, 56221 
(2001). 

[40] L. E. Reichl, G. B. Akguc Found, of Phys. 31, 243 (2001). 

[41] Bryce S. DeWitt, Rev. Mod. Phy. 29 377 (1957) 

[42] William H. Press, Saul A. Teukolsky, William T. Vet- 



tcrling, Brian P. Flannery, Numerical Recipes in For- 
tran.: The Art of Scientific Computing Second Edition. 
Cambridge University Press 1992 

[43] M.C. Gutzwiller, Chaos in classical and quantum me- 
chanics, Springer 1990. 

[44] Leonid I. Glazman and Michael Pustilnik, Nanophysics: 
Coherence and Transport eds. H. Bouchiat et al. (Else- 
vier, 2005), pp. 427-478 



